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Abstract 

A  calculation  of  the  ionization  cross  section  and  the  excitation  cross 
section  to  the  2S  and  2P  states  of  hydrogen  by  collision  with  electrons  is 
carried  throxigh,  using  an  impulse  approximation.  The  resxilts  are  then 
conrpared  to  the  experimental  results  when  available  and  to  the  various 
theoretical  calculations.  The  cross  sections  obtained  compare  favorably 
with  experiment  except  for  the  excitation  to  the  2S  state.  The  calcu- 
lations carried  out  by  these  methods  are  considerably  simpler  than  the 
usual  approximation  methods,  and  the  agreement  of  the  results  with  ex- 
periment indicates  that  it  is  a  simpler  and  better  method  to  use  for 
the  computation  of  cross  sections  in  the  scattering  of  electrons  by  atoms. 
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1.   Introduction 

The  problem  of  the  scattering  of  electrons  by  atoms  is  of  theoretical 
interest  since  it  serves  as  a  model  for  other  many-body  scattering  pro- 
cesses and  since,  in  this  case,  the  law  of  force  between  the  particles 
involved  is  known.  Ihere  is  also  a  practical  motivation  in  this  problem. 
Whenever  one  considers  the  passage  of  an>  electron  through  a  cloud  chamber, 
the  propagation  of  electromagnetic  reidiation  through  the  ionosphere  or 
gaseous  discharges,  the  scattering  of  electrons  by  atoms  plays  a  pre- 
dominant role.  Inelastic  processes  play  an   important  role  whenever,  in 
a  physical  process,  an  electron  transfers  part  of  its  energy  to  the 
internal  degrees  of  freedom  of  the  scattered  system.  The  pvirpose  of 
this  paper  is  to  apply  an  approximation  method,  hitherto  unused  in  this 
field,  to  the  analysis  of  this  problem. 

We  shall  be  dealing  in  this  paper  with  the  ionization  of  atomic 
hydrogen  and  the  excitation  of  a  hydrogen  atom  from  its  ground  state  to  the 
n  =  2  level  due  to  a  collision  with  an  electron.  Our  interest  in  hydrogen 
is  two-fold.  Atomic  hydrogen  is  one  of  the  basic  constituents  of  the 
universe,  and  thus  the  problem  has  an  intrinsic  interest  of  its  own. 
Beyond  this,  the  scattering  of  electrons  by  hydrogen  can  be  used  as  a 
model  for  approximation  techniques  which  caxi  possibly  be  generalized  to 
deal  with  more  complex  situations.  In  recent  years,  experimental  techniques 
have  improved  so  much  that  one  is  able  to  cariy  out  experiments  in- 
volving atomic  hydrogen  and  analyze  the  results  to  separate  the  various 
inelastic  processes.  For  exaiqple^Fite  in  the  United  States  and  Boyd  and 
Green  in  Bigland  have  been  able  to  investigate  the  ionization  and  ex- 
citation of  hydrogen  by  electron  impact.  In  addition  to  yielding  data 
with  regard  to  the  elementary  processes  involved,  the  results  serve  as 


a  test  for  the  various  approximation  methods  used  for  the  analysis  of  the 
inelastic  scattering  of  electrons  by  hydrogen  atoms. 

The   approximations  that  have  been  used  to  date  are  the  Bom  approxi- 
mation^ the  distorted  wave  method  and  some  variational  techniques.  The 
problem  Involved  in  an  approximation  method  is  to  find  a  reasonable  des- 
cription of  the  wave  function  of  the  electron-hydrogen  system.  In  the 
Bom  approximation,  one  assumes  that  the  interaction  between  the  electron 
and  the  hydrogen  atom  is  so  small  that  one  can  choose  for  the  total  wave 
function  the  product  of  a  plane  wave  for  the  incident  electron  and  a  wave 
function  describing  the  xmperturbed  atomic  system.  In  the  distorted  wave 
approximation,  one  tries  to  make  an  estimate  of  the  distortion  of  the  in- 
cident wave  by  the  atomic  system.  In  using  variational  methods,  one  cal- 
culates some  expression  describing  the  scattering  process  which  is  station- 
ary with  respect  to  the  choice  of  the  xuiknown  wave  function  of  the  system: 

that  is,  one  hopes  that  an  error  of  order  e  in  the  wave  function  will  re- 

2 
suit  in  an  error  of  order  €  in  the  description  of  the  scattering  process. 

A  coniparison  of  the  experimental  results  with  the  theoretical  results 
shows  that  the  Born  approximation  tends  to  overestimate  by  a  considerable 
amoxmt  the  value  of  the  cross  section  at  low  energies.  In  general,  the 
maximum  of  the  calcvilated  cross  section  occurs  at  a  value  of  the  energy  less 
than  the  value  at  which  the  actual  experimental  maximum  occurs.  Even  the 
second  Bom  approximation  gives  a  value  which  is  still  too  Isirge  for  the 
cross  section,  although  it  is  closer. 

In  this  paper,  we  have  calculated  the  inelastic  cross  section  to  the 
n  =  2  level  and  the  ionization  cross  section  of  hydrogen  for  collisions 
vlth  electrons  using  the  'Impulse  Approximation*  developed  by  G.  Chew '-J. 
Ihe  basic  Idea  of  this  approximation  is  to  consider  the  target  particles 
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as  frozen  dixring  the  collision  -  that  is,  to  disregard  the  fact  that  the 
atomic  particles  are  bound  rather  than  free.  If  the  cross  section  for  a 
collision  between  two  particles  is  known,  as  it  is  in  this  case,  then  one 
can  calculate  the  cross  section  due  to  the  two  two-body  collisions.  To 
take  into  account  the  fact  that  the  two  target  particles  are  bound,  one 
averages  over  the  charge  distribution  of  the  bovmd  particles  in  the  initial 
and  final  states.  In  order  to  account  for  the  assumption  that  the  boimd 
particles  do  not  interact,  one  replaces  the  charge  distribution  with  an 
equivalent  momentum  distribution. 

The   validity  of  this  approximation  can  be  seen,  from  physical 
considerations,  to  involve  the  transit  time  of  the  incident  particle 
across  the  target  system  and  the  period  of  the  bound  particle.  To 
estimate  the  circumstances  under  which  this  approximation  is  valid 
for  the  scattering  of  electrons  by  hydrogen  atoms,  6ne  proceeds  in  the 
foliar ing  way.  The  radius  of  the  hydrogen  atom  in  the  ground  state  is 
a  .  If  the  orbital  velocity  of  the  bound  electron  is  v  ,  then  the 
period  of  the  electron  is 

If  V-  is  the  velocity  of  the  incident  electron,  then  the  time  of  transit 
is 

t-  =  2a  /v,  . 
1      o'  1 

The  validity  of  the  inrpulse  approximation  depends  on 

*1  «  *2 
or 


Vg/Vj^      «      Jt. 


-  ll . 

At  threshold,  v.  =  v  and  thus,  as  one  approaches  threshold,  the  results 
nay  not  he  too  reliable.  As  v^  increases,  v-/vj^  decreases  so  that  for 
an  incident  energy  of  100  e.v.  v./v.  =  I/5.  Hence,  as  the  Incident  energy 
increases  the  approximaticn  should  yield  resvilts  that  cortrpare  more  favor- 
ably with  experiment. 

This  paper  is  divided  in  the  following  way:  a  review  of  the  various 
approximation  methods  is  taken  up  in  Section  2.  In  Section  5,  we  derive 
the  fundamental  fonm.ilas  of  the  Impiilse  Approximation  and  state  these 
for  ionization  and  excitation.  Section  k   deals  with  the  ionization  cross 
section  while  Sections  5  and  6  deal  respectively  with  the  cross  sections 
for  the  2S  and  2P  states.  Section  7  summarizes  the  results.  The 
evaluation  of  the  various  integrals  occurring  in  the  body  of  the  work 
is  done  in  Appendices  A,  B  eind  C. 

2.   Review  of  former  work 

The  problem  of  the  scattering  of  an  electron  by  a  hydrogen  atom  can 
be  found  as  a  solution  of  the  following  Schr'ddinger  equation  with  appro- 
priate boTUidary  conditions: 


(1)  -!  -  a  (^  ^  '^2>  -  fj  -  fj  -  4  ^  5  (^i'^2'  =  ^  5  ('r^a' 


In  the  above  equation,  the  center  of  mass  of  the  system  is  located 
at  the  proton  which  is  assmned  to  have  an  infinite  mass;  r^  and  r^  are  the 
coordinates  of  the  incident  and  bound  electrons,  respectively,  and  E  Is  the 
total  energy  of  the  system.  The  equation  is  not  separable  because  of  the 
term  e  /r,p-  We  shall  now  discuss  the  boundary  conditions  to  be  imposed 
on  the  solution  of  (l)  and  the  information  which  is  needed  from  the  solution 
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Later  ve  shall  discuss  the  various  approximations  that  have  been  used. 
We  desire  a  solution  of  (l)  which  has  the  asymptotic  form 

Ik  "^  ^n^l 

r.— >oo  J. 

where 

and  ^(Pp)  are  the  hydrogenic  wave  functions.  The   svmimatlon  includes 
integration  over  the  continu\im  states,  f  (9-,0.)  is  the  scattered 
amplitude  for  excitation  to  the  nth  state.  It  is  related  to  the  differ- 
ential cross  section  for  direct  scattering  by  the  relation 

(5)  <^nd.^^>^^^  "^   =  k  K^^'^^  1^  *"• 

o 

Instead  of  taking  the  limit  as  r^  — >  od  ,  we  can  allow  r^  to  approach 
00 .  Ihis  win  describe  the  process  wherein  the  Incident  electron  ex- 
changes roles  with  the  bound  electron.  The   asymptotic  form  of  the  wave 
function  in  this  case  is 

W  1  (^^.r^)  ^    r  g„(92^02>V^l)  ^-T— 

rp->  00  n  2 

where  zA^pt^o)   is  the  exchange  scattered  amplitude  and  is  related  to  the 
exchange  differential  cross  section  by 

(5)  ,      «;je,0)  dm  =gS  \g^iQj)\\m     . 

o 
Ihe  dlfferentieil  scattering  cross  section  for  the  process  is  a  ccoibinatlon 


-  6  - 
of  direct  and  exchange  scattered  amplitudes .  One  must  also  take  into 
account  the  spins  of  the  two  electrons.  This  leads  to 


o 


Let  us  now  derive  the  general  relations  upon  which  most  of  these 
approximations  are  based.  To  this  end,  we  expand  ^  (r^,r  )  in  terms 
of  the  complete  set  of  hydrogenic  wave  functions  ^  (r-) 

(7)  lir^.r^)     =  n^ir^Kir^)' 

As  usual,  the  sximmation  sign  includes  an  integral  over  the  continuvmi 
states  of  hydrogen.  Substitution  of  (7)  into  (l)  leads  to 


(8)         EI„(?,)^  ^  -  ^^^F„(?,) .  ^J  ^  .  U|  (?^,?^) 


where  we  have  used  the  Schrodinger  equation  for  the  hydrogen  atom  satisfied 
^y  5  (^p).  We  multiply  both  sides  of  (8)  hy  ^  (r^)  and  integrate  over  r-. 
In  making  use  of  the  orthogonality  of  the  hydrogenic  wave  functions,  we 
obtain: 

(9)      iV  ^nl^  V"l)  =  21  V^l^V'l^   ("  =  0,1,2,...) 


where 


(10)      U      ^^ 


nm     4-2 


t' 


^aO^^M^-HlJ?,) 


,-i   -=>> 


For  exchange  scattering,  we  expand  ^  (r^,rp)  in  terms  of  the  complete  set 
of  hydrogenic  wave  functions  $  (r. ) 


(11)  l^tr^,?^)     =     E0,(?2>l„(?i)- 
Application  of  the  same  aiethod  used  to  derive  (9)  leads  to 

(12)  (A^.  ^l)  G(?,)  =  ^  j  ^1  ^n<^i%  -  ^)  ^  '^'^2'   • 

Each  of  equations  (9)  and  (12)  consists  of  an  infinite  number  of  coupled 
integro-differential  equations  and  cannot  be  solved  exactly.  They  are 
solved  approximately  by  using  various  schemes.  We  shall  look  into  these 
briefly. 

A.   Born  approximation 

Let  us  first  consider  (9).  The  Born  approximation  is  obtained  by 
assimilng  that  all  the  terms  on  the  right-hand  side  of  (9)  are  negligible 
except  for  the  term  associated  vith  the  incident  wave.  This  is  equi- 
valent to  assuming  that  the  pertvirbation  due  to  the  interaction  between 
the  two  electrons  is  small.  Equations  (9)  take  the  form: 

(15)  (A,  +  k^)  F  (rj  =  U  e  °  ^   . 

^   '  ^1    n'  n^  1'    no 

2 

To  solve  this  equation,  we  use  Green's  function  for  the  operator  A.  +  k 

satisfying  the  appropriate  boundary  conditions.  This  yields 


,r^  -»•    ik  r,  -  r' 
(IM     F„(r,)  -.^i   U^J.,)e  °  1  ^  ^  d?i  . 


=■1  -  n' 


If  we  now  teike  the  limit  as  r^  — >  00 ,  we  obtain 

^n^l 
(15)  F^(r,)       -_^  ^   fjBj)  ^ 


with 
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(16)         f  je,0)  =  -  J^  f  U^Jr,)e'^^°"^°^"^  d?^ 


For  Ionization,  (l6)  takes  the  form: 


(17)  v»,^)  =  - 1;  ^  J  a?,^^  /"'-^^  i;  (;,)(i  -  ^je"""  •  ^1  j„(?3) 

In  equation  (17)  \     (r  )  is  the  vave  function  of  an  electron  of  momentum 

In  the  continuum  field  of  the  proton.  Mott  and  Ma8sey'--J  use  the  expression 

where 

i   =  r(l+co8  0  )      n  =  l//<  a  . 

cos  0  =  COS  e^  COS  e^  +  sin  0^  sin  9  co8(e  -  0  ) 

and  J^(2  /iK.  |u  )  is  the  Bessel  function  of  zero  order  and  p(l-in)  is 
the  gamma  function.  If  we  substitute  (iS)  Into  (l?)  and  integrate  over 
r^  and  r^,   we  obtain  for  the  differential  scattering  cross- section 


a  kK^ 
o 


(19)  exp{-2ft<„-l(2^^/  [k2.^V])1 

"  {.^  .  a.^,A  ^^).(K=-  .^)^j  '     ""^ 

where  integration  over  an  angles  of  ejection  of  the  bound  electron  has  been 
carried  out.  In  (19)  1^  is  the  wave  number  of  the  ejected  electron 

M  =  Z=l/a^ 

K^  =  (k^  +  »\.^  -  2k  V^^cos  ©) . 

This  expression  was  integrated  numerically  by  W.  Fite  ^^ .     The   ciorve  obtained 


is  plotted  in  Figiire  1.  Note  that  at  high  energies,  the  value  of  the  total 
cross  section  is  too  high. 

For  the  excitation  to  the  state  n  from  the  state  m  by  electron  impact, 
Mott  and  Massey  L  -I  give  the  expression  for  direct  scattering 


(20) 


W^^*^ 


aV  k5  ''-^^^' 


where 


(21) 


mn 


e     J  5  dr  . 


For  excitation  from  the  ground  state  to  the  state  (nj2m)  they  give: 


c  ,„  =  22^*V+^2/+i)^/2(t+i): 

n  JLID. 


(n-i-l)i 


(iH-^)! 


1/2 


(Ka^) 


(22a) 


X 


f(n-l)^44C^} 
{(nfl)2+4^2] 


l/2(n-/-5) 
l/2(n+X+3) 


{r^l)\{n-lf-^kAc^\_^M 


-   2n'((n-l)^44C^ 


^^^',.x,=.c4^^V„-.,(.. 


+  (n-l)J(nfl)2+4^2l  c^+2^_^(x) 


where 


(22b)     X  =  (n^-1+45^) 


l(ttfl)2+4?2U(^.^)2^^2l 


-1/2 


(22c)     ?  =  |Ka^ 


and  C  (x)  are  the  Gegenbauer  polynomials.  The  ctirves  obtained  for  the 
s 

total  cross  section  for  excitation  to  the  2S  and  2P  from  the  ground 
state  are  plotted  respectively  in  Figures  2  and  3.  Again  the  values  are 
too  high. 
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Rothenstein^-  J  has  carried  out  the  second  Bom  approximation  for  ex- 
citation to  the  2P  state.  Although  the  values  are  better  at  intermediate 
energies,  they  are  still  not  in  good  agreement  with  experiment.  Ihese 
results  aure  plotted  in  Figure  5« 

B.   Distorted  wave  approximation 

In  this  method,  one  assumes  that  all  the  non-diagonal  elements  of 
the  matrix  U   given  in  (9)  are  small  except  for  the  ones  associated  with 
the  incident  wave.  No  assimrptions  are  made  regarding  the  diagonal  elements. 
Tbaxs   (8)  reduces  to 

(25a)         (^  +  ^^0-  W"l>  =  0 


(25b)         (^^J^n-  W^l)  =  "noV^l^  '  ^(^1^ 


where  H(r- )  is  a  known  function.  If  the  solution  of  the  homogeneous 
e  quatlon 

(25c)         (^  +  ^^n-  W^l^  =  ° 
is  \r{r^,a-<S)),   then  the  solution  of  (25b)  for  large  r^  is 

ik^r^ 


The  scattered  amplitude  will  then  be 

(25)   ,       f^(e,0)  =  -  J;j  I  d?i  i|^(ri,«.  ©•)  U^o(ri)Fjrp 

irtiere  ^^(r^)  is  the  regular  solution  of  (25a)  and  U    (r)  is  given  by  (lO) 


-  11  - 

Massey'-  -'  points  out  that  in  order  to  include  exchange  in  this  approx- 
imation, one  can  use  the  close  coupling  treatment.  In  this  case,  instead 

-_   -1.   -3» 

of  using  (6)  as  an  expansion  for  the  wave  function  J  (r^^r.),  we  use  a 
different  form.  Since  we  are  considering  two  electrons,  ^  (r^,r2)  is 
either  symmetrical  or  antisymraetrical  depending  on  whether  the  spins  of 

the  electrons  are  parallel  or  antiparallel .  We  assimie,  then,  that  to  a 

_  _i.  ^ 
good  approximation,  J  (r-,rp)  is: 

(26)      1  (r^,^,)  =  F*(r^)^^(r^)  .  Fp^)^^^^)  .   F*(^,)1J?,)  ±  F^Cr^)! J?^)  . 

The  calculation  has  been  applied  to  the  excitation  from  the  IS  to  the  2S- 
state  of  hydrogen.  Massey'-  -i  shows  that  in  this  case,  one  is  led  to  the 
following  expression  for  the  scattered  amplitude 

fti^J)  =  -  fepi(-r«-®i)V^i)^o(^i)n 

(27) 

*fe  j'^i4^i(-r-®i)Vv-2)^o(^2) 


where 


v^i)  =  ^/o^v5i(^2)  ?^-B 


(28a)         U,,(rJ  =  2  |  ^:(r,}$,  (r  J^  ^  -  ^^  dr^ 


(28b)         ho^r^'-2^   =  h^^A^^2^H   "  ^  ~  '"^o 


and  •v/   and  ^i^   are  distorted  waves.  In  order  to  proceed  further,  Massey*-  -' 
o       1 

carries  out  a  spherical  harmonic  analysis  of  v/  and   v^  and  points  out 

that  only  the  zero  order  cross  section  has  been  calculated.   In  this  case, 

the  equation  satisfied  by  */    is 

o,o 
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with  a  similar  equation  for  J  ^  f.-     Here 


K°^(r,--)  =  4nrr'^(r)l^(r'K-27Q(r,r')  +  k^  -  2E 


where 


7  (r.r ' )  =  — r      r'  >  r 
'  o  r 


-  .      r'  <  r 

r 


The  functions  -^  and  ^/ ,  ^  cannot  be  obtained  directly  and  niust  be 

o,o      1,0 

obtained  by  assuming  a  trial  function  which  asymptotically  behaves  as 

»±  1  _... 


.y1*  Aj  A*  -  sin(k  r  +  ■)?*) 
^  o,o     or    ^  o     ^  o' 

^.^*  ,  ^  B*  i  sin(kTr  +  (,*) 
0,1      or    ^1    ^o 


where  f?     and  ^  are  the  zeroth  order  phase  shifts.  The  resiilts  of  these 

calculations  show  that  a  maximum  occurs  at  about  15-5  e.v.  However,  the 

2 
magnitude  of  the  calculated  cross  section  as  given  by  Massey,  O.Tl^fla  , 

-2  2 
is  about  ten  times  larger  than  the  experimental  value,  6  x  10  jta  . 


3.   Impulse  approximation;  fvuidamental  theory 

We  now  derive  the  various  formulas  we  need  for  our  calcxilations .  Ihe 
Schr5dinger  equation  for  the  problem 


(1)  H  1  (r^,r2)  =  E  f  (r^^Tg) 


where  r.  and  r^   are  the  position  vector  of  the  incident  and  boimd  electrons 
respectively  and  E  is  the  total  energy  of  the  system  given  by: 


(2)  E  =  ^-  |e| 

_^ 
where  k^  is  the  wave  number  of  the  incident  electron  and  e  the  binding 

energy  of  the  groimd  state  of  hydrogen.  H  is  the  Hamiltonian  of  the 

system 

t2     x2      2    2    2 
(5,  H  =-|jA,-fcA..^-^.^ 


2m  2   r^  .r^   r^^ 


We  desire  a  solution  of  the  form 


Here  ^  .  (r^^Tp)  represents  the  initial  state  of  the  system  and  is  given  by 

ik  'T 
(5)  li(V^2^  =  (2rt)-5/2e  ^'   ^I^ir^) 


ik^-r^ 
where  e      represents  electron  (l)  as  an  incident  plane  wave  of  momentum 

k^  and  $.   (r-)  is  the  wave  function  describing  the  groiind  state  of  hydrogen. 

5  (r  ,T   )   represents  the  final  state  of  the  system  and  contains  various  kinds 

of  outgoing  waves  and  is  given  by 


(6) 


lj^^l'^2^   =  J  G(?^,?2'^|,?p  V  f   (?^,rpd?^  d?^ 


'^  -^  -^  -^  .  I 

where  G(r^,rp;r',r' )  is  the  Greens  function  for  the  equation 

(7)     (H^-^)  K'^.r^)  =  0 
with 

V2      i.2      2 

'«  =0  =-1:  ^-s^2-^    ■ 

A  Green's  function^  -•  for  our  problem  is 

ik' |r  -r^l 

(9)        G(?,,;3,?i,?p  =  -  ^  E  J,(?2)i(?|)  ^  !     .,, 

2jm   n  I  r  -  r  •  I 

where  the  svmi  includes  an  integral  over  the  continuum  states  and 


■Ik- 


(10)       -^l   '   p(=-\)- 


For  large  r^ ,  we  then  have,  upon  substitution  of  (9)  into  (6), 
Is  <?1'^2>  r2^>  <.  -  ^  r  5.(?a) 


n  1 

e 


(^)  .  -ik'.?: 


d?{  <i?^  5^(rpe   '^  ^/  f  ( r^,rp  . 


For  exchange,  we  tpke  the  limit  as  r^  approaches  oo,  and  this  leads  to: 


ik'r. 


(12) 


/  -ik  -r 

J  d?^  d?^|  *(?pe   °  ^  $  (?{,rp. 


In  these  equations,  ^  represents  the  wave  function  of  the  electron  in 
the  field  of  the  proton  in  the  state  n.  We  assrune  the  following  approxi- 
mate form  for  ^  (r-,r_) 

(13)      $  (r^^r^)  ->  I3  (^i^r^)  =  f  di^  gjk^)  f  t^^t^irj^)- 


Here  g  (k^)  is  the  Fourier  transform  of  the  hydrogen  atom  ground  state 
wave  function.  $  ^  ^  (r-,r„)  is  a  wave  fvinction  representing  the  Inter- 
action  of  an  electron  of  momentum  k^  by  a  free  electron  of  momentimi  kp. 
What  we  have  done  in  (I5)  is  to  replace  the  bound  electron  by  a  free 
electron  whose  momentimi  distribution  is  equivalent  to  the  ground  state 
wave  fiinction  and  consider  its  interaction  with  the  incident  electron. 
liie   substitution  of  (15)  into  (12),  gives 
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(1^)  ^,_^, 


Por  exchange,  we  take  the  limit  as  r_  approaches  od 


2 


(15)  ff    ^    ^      -i^^-r^        _.^ 


J 


dr^dr-e   "  "  InC-^)^  1  ^^^^(^i^^P 


We  can  represent  the  scattered  amplitude  by  an  R-matrix.   In  our  case, 
this  is 


-ik'-r,    .  ^  ^ 


(16)      (k',n  lR^2l^i>^2^  =  -^  Jj  ^1<^^2  ^   "  '^n(^2)^  1  t^ij-'^l'^2^ ' 

The  R-matrix  represents  the  matrix  elements  of  the  perturbing  potential 
between  the  initial  and  final  states  of  the  process  and  is  related  to  the 
scattering  cross  section. 

In  order  to  derive  the  relation  between  the  R-matrix  and  the  scattering 
cross  section,  we  use  the  time -dependent  scattering  theory  developed  by 
Gell-Mann  and  Goldberger  L"j ,  Our  system  is  described  by  the  following 
equation  (^  =1) 

(IT)  i-^l(t)  =  H$(t) 

St 

where  H  is  the  Hamiltonian  and  ^(t)  the  wave  function  of  the  system  composed 
of  an  electron  and  a  hydrogen  atom.  We  separate  H  into  tv^o  terms  by 

(18)  H  =  \  +  \ 
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where  H^  is  the  Hamiltonian  describing  the  system  when  the  interaction  between 
the  incident  electron  and  the  hydrogen  atom  is  cut-off  and  H^  represents 
the  interaction  terms.  Ihus, 

2         2 
(I8b)      H  =  -  f-  +  f-    . 
'^      1    12 

The  non- interacting  state  can  then  be  represented  by  the  equation 

^(t) 
(l8c)      i  -5^—  =  H^l(t)  . 

The  purpose  of  ^(t)  is  to  characterize  $('*')• 

We  are  interested  in  transitions  from  the  inital  state  ^.  to  the  final 
state  5..   If  the  system  was  in  the  state  ^.  at  an  infinite  time  in  the  past, 
we  can  assume  that  J.  can  be  described  by 


(19)      li(t)  =  e-^i^(T-^^$i(T) 


Equation  (I9)  implies  that  the  incident  wave  is  released  suddenly  into  the 

system.  Actually,  the  beam  is  prepared  in  some  way  and  is  fed  into  the 

.system  over  a  period  of  time.  To  account  for  this  fact,  instead  of  using 

(19),  we  use  a  statistical  average  of  waves  of  the  form  (I9)  over  a  period 

of  time  T  and  take  the  limit  as  T  ->  oo.  This  leads  to 
o  o 


(20)      l.(t)  =  liin     ^  [  dTe^^C^-*^^  (T) 
^      T  ->  00  o  >„ 


^o 


Writing  e  =  1/T^  , 
.(e) 


(21)     1^  (t)  =  e  j  dT  e'^  e^^C^-*^  ^^(T)  , 

-  00 
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remembering  that,  in  all  evaluation  of  the  cross  sections,  ve  take  the  limit 

T  '  5 

as  e  ->  0.  The   ^  s  are  normalized  in  a  cube  of  side  L  and  volvune  L  .  Since 

T  =  L/v  where  v  is  the  group  velocity,  it  follows  that  the  limit  €->  0 
o 

and  L  ->  oo  must  be  taken  simultaneously.  Furthermore,  as  e->  0  and  L  ->  oo, 
e"-4."^  ->  0. 

We  must  now  derive  the  rate  at  which  transitions  ai-e  taking  place  from 
the  state  i  to  the  state  J .  The  probability  for  transition  from  a  state  i 
to  a  state  j  is: 

where 

(22a)    f^j(t)  =  <lj(t)||^(t)> 

(22b)    Nj  =  <$j(t)|$^(t)>  . 

Let 

-iE  t 
\{t)   =  e   i  $. 

then  (l8c)  gives 

(23)    (H^  -  E^)  1^  =  0  . 


Equation  (2l)  then  becomes; 


-  CD 


(..)  ilM  -  .  (dT  .--e%  ^  ■  ^'  1,  =  e:.-^  ^-^^  1, 


But  ^.  is  an  eigenf unction  of  E.  and,  if  we  consider  H  operating  on  J  ,  we  have 
e\   =  (H^  +  H^)  $. 
=  (E^  +  H^)  l^ 


or 


(25)     (H  -  E^)  li  =  Hg^.  . 
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•Hien  at  t  =  0, 

or 

(26)    lf^>(0)  -  \  *  5-^  Hj^ 

using  (25).  Instead  of  (26),  we  have,  "by  solving  for  "§. 
*i   %    E.-H^+ie  % 


or 

(27) 


i 


(28) 


^   J 

vhere 

(29)       Rij(e)  =  <lj|H2ll^(0)  >  . 

If  we  expand  (27)  in  a  power  series  of  ^.,  we  see  that  the  term  „  _„     . 
acts  ohly  on  V^  and  not  on  J  directly.  Itus,  R   is  a  well-behaved  function 
of  the  energy.  Further,  since  R..  is  proportional  to  L  ,  due  to  the 
normalization  of  the  \  s,   it  follows  that  as  e->  0  and  L  ->  00, 

(50)     R.j(e)L5->R.^. 

We  mentioned  above  that  we  were  interested  in  the  rate  of  transitions,  that 

•^^ij 
is,   J."  .  To  calculate  this,  we  need  the  time  derivatives  of  f. ..  We  write 

i(E  -H)t  /  X 
(51a)     fj^j(t)  =<lj|e   ^    lli'^(0)>. 
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Then 


(51b) 


-^ =   i  <$j|Ej   -   H|f  ^^0)  > 


dt 


(O. 


KiJHplli'^CO)  > 


i  R^j(e)    . 


Since  from  equation  (22), 
dW 


11 


dt 


.=r  ^  t^  '^«<^' 


t=o 


we  must  calculate  —  |f..(t)|  .  We  have 


—  If  (t)l^ 
dt'^ij^^^l 


t=o 


=  fij(O)  flj(o)  .  f,j(0)  f*.(0) 


where  the  dot  denotes  differentiation  with  respect  to  time. 
Using  (28)  and  (51),  we  obtain 


(52)      t^ifij(t)r 


=  26^jImR^j(c)  + 

Jt-O  v-j  ^^ 


iJ— iJ^^'  ■  (E  -Ej2^e2  '  iJ 


R..(or 


We 


must  now  deal  with  N. .  The  ^  s  form  a  complete  set  of  states.  Hence,  if 


we  compute  YZ   |f^*('t)|  from  equation  (51a),  we  obtain 


O^) 


z:if,,(t)i^  =  N^ 


Furthermore,  N.  is  a  normalization  constant  and  hence  independent  of  time. 
It  follows  that  if  we  take  the  time  derivative  of  (55)  at  t  =  0,  we  have 


(5^^) 
Thus 


2  ImR  (€)  +Y1 

'^'^  j   (E.-E,)'+€ 


^H:2  iRij(oi^  =  o 


i  J' 


N,  =  1  +  -  I  R.,(e) 


by  using  (28)  and  (5^).  Now  as  €  ->  0  and  L  ->  od  we  see  that  N  ->  1. 


Tlius,  the  transition  rate  is 


(55) 


W. ,(0)  =  26. .1  R.  .(c)  + 


2e 


J  1 


R..(or  . 
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Ihe  scattering  cross  section  is  equal  to  the  transition  rate  divided  by  the 

-5 
initial  flxuc.   In  our  case  the  initial  fliix  is  vL   where  v  is  the  velocity 


of  the  incoming  electron.  Hence, 

"■^        €->0  (E-E  )^+^   ^-^ 
L->OD 

As  e->  0,  the  factor  2g/(E  -E^)^+G^  ->  2n6(E  -E^)  ,   and  by  using  (50)  we  obtain 

This  is  the  relation  between  the  R-matrlx  and  the  cross  section  that  we  sought. 

Let  us  go  back  to  equation  (l6).  The  potential  V  consists  of  two  terms 
2      2 

^1     ^12 

Substituting  into  (l6),we  obtain       _^ 

/  /■       -  ik '  I" 

If  we  were  to  use  the  exact  wave  functions  for  the  bound  electron  in  the  final 
state,  i.e.,  Coulombic  wave  functions,  the  integration  over  r-  would  yield 
zero  when  we  compute  the  matrix  element  involving  l/r^  since  the  wave  functions 
in  the  initial  and  final  states  are  orthogonal.  Tius,  it  would  seem  that  the 
wave  originating  at  the  proton  does  not  contribute  to  the  inelastic  scattering; 
the  only  effect  of  the  proton,  which  we  assume  to  have  Infinite  mass,  is  to 
deflect  the  incident  electron  into  a  new  direction.  Thus,  in  inelastic 
processes,  the  wave  originating  at  the  proton  does  not  contribute  to  the 
scattering  cross  sections . 

One  further  simplification  will  be  made.  Instead  of  using  hypergecanetric 
functions  to  describe  the  electrons  in  the  continuum)  ,  we  will  use  plane  waves 
for  ionization.  The  Integrals  are  intractable,  otherwise. 

One  immediate  effect  of  this  last  siBrpllfication  ajta  be  seen  very 

ik  -r 

quickly.  Looking  back  at  (l6)  and  putting  l„(r  )  =  e  ^  2  we  see  ths.t  the 
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R-matrix  in  the  case  of  direct  scattering  is 

-  ik '  •  r  -  ik '  •  r 

while,  for  exchange, 

-  ik^  T^  -  ikA  •  r, 

■  vl- 

1"2 


(4,k' IrJ^^)  |k^,k2)  =  -  i  jj  ^^^^e'  "^^   '2   ^^2  •  ^1  ^  ^  g.  (-^,?^) 


Ihese  two  expressions  are  exactly  the  same  except  that  we  have  interchanged 
r,  and  r  or,  by  the  same  token,  k.  and  k  .  This  will  then  have  no  effect 
on  the  final  value  of  the  integrals.  The  two  electrons  are  being  treated 
identically  axid  we  cannot  distinguish  between  singlet  and  triplet  scattering. 
The  main  difference  between  the  correct  wave  function  and  a  plane  wave  is  the 
fact  that  in  the  latter,  the  singularity  at  the  origin,  in  the  wave  function, 
is  absent.  The  effect  of  this  approximation  has  not  been  estimated.  Since 
our  restilts  compared  favorably  with  the  experimental  data  available,  we  rely 
upon  these  for  our  justification  in  using  this  approximation. 

We  factor  J  v  v  ^^i^o)  ^^^°   *^®  product  of  two  functions,  one 

k^kg  1  d 

representing  the  motion  of  the  center  of  mass  sind  the  other,  the  relative 
motion  *ith  respect  to  the  center  of  mass: 


Ihus,  (l6)  becomes 

Substituting  for  ^(r-)  its  Fourier  transform,  we  get 


ICk^-k^) 


-  22  - 


(57)    (4,n,lR3^2l^l'^2^ 
where 


=  e 


-|dk''g^(k|)(4,k''iRi2l^i,t 


(58)    (k',kj|R^2l^r^2>  ^  "  (t?  Ji  ~tr  "^   pi'V^^2  -V  ii(V'^2)-(V^2) 


^(k^-?^) 


(Vrg). 


We  now  caxry  out  the  transformation 


r  -  r  =  r 
1    2 


2(^1  +  r  )  =  p  . 


i'^2  l^2l^r^2^  =  ■  "^^  J  <i?exp,  -i(4+k^'  -k^-kg)-?.  I 


nxis  leads  to 

(39) 

where 


-  e      Vr) 


=  5(k^  +  k^'  -  k^  -  k^)(k'  Irjk) 


k  =  I  (k^-k^) 
k-  =  I  (k^-k^'l 


and 

(ho)       (k* 


Hie  R-matrix  describing  the  inelastic  collision  will  then  be 


{kl)  (k|,n|R|k^,0)  =  jj   dk''dk^"g*(k^')(k^,%  I  Ri2i^i'%')eo(%')  • 

Substituting  (59)  into  (i+l)  and  carrying  out  the  k"  integration  leads  to 

{k2)  (k^n|Rjk^,0)  .=  e^j  dk^'  ^(k^O(k'  |r|k)g^(k^+k^'  -  k^) . 

For  ionization,  the  R-matrix  is  obtained  by  letting  n  correspond  to  k'. 
g^(kp')  will  then  be  the  Fourier  transform  of  a  plane  wave.  Tbus 
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^2 
For  excitation  to  the  2P  and  2S  states,  equation  (42)  is  applied  directly. 
For  exchange,  we  see  from  (12)  that  in  the  case  of  ionization,  the  only  effect 
is  to  interchange  r^  and  r  in  the  final  state.  Since  we  are  using  plane 
waves,  there  is  no  distinction  between  the  two  scattering  processes.  For  the 
2S  and  2P  states,  we  have 

(44)   (k^  ,  n!R^|k3^,0)  =  -  ie^J  dk2go(k2)g^(k^  +  k^  -  k^)  |  f^  e^^-\it). 

Before  we  proceed  to  compute  the  cross  sections  we  consider  the  r-matrix. 

y  (4( 

.[10] 


This  is  given  by  (4o).  H  ^r)  represents  the  interaction  of  two  free  electrons 


and  is  given  by' 

(45)  E^{t)   =  N(k)e^*^  F(-  i/2a^k;l;ikr  -  ik-r) 

Wliere  F  is  the  hyper  geometric  of  its     argument  and 

(46)  N(k)  =  ^  fd  +  i/2a^k)e-  «/^^o^  > 

where    P(l+i/2a  k)   is  the  garama  function 
Thus   (40)  becomes 

(47)  (k'|r|k)  =  -  iN(k)l 
with 

(48)  I  =    I  f^  gi(k-k')-r  p^_  i/2a^k;l;ikr  -   ik-?^    . 

We  rewrite  this  integral  as 

-   -V.     -k  -fc 

(49)  I  =  J  F"  ^^"^'^  ^^"  i°>l^^^  -  ^•^• 

NordsieckL  J  has  evaluated  this  integral  for  the  case  in  which  an  exponential 
factor  of  the  form  e~   occurs  in  the  integrand.  We  have  evaluated  this 

integral  in  Appendix  A  and  the  result  is 

I    2   ^-^  f  in 

(50)  I  =  ^  (^--^)   . 

<1     q 


2k  - 


We  will  now  proceed  to  apply  the  above  formulas  to  the  cases  we  have  under 
consideration . 

k.        Ionization  cross  section 

The  R-matrix  for  ionization  is  given  in  equation  (jAj)  which  states 
that 
(1)    (4,k^iRi^,j,l%0)  =  e^  I  dk^  (k'  |r|k)g*g',  (k^'  )g^(k5^  +  k^'  -  kj^) 

with  (k'|r|k)  given  by  (5  AT)  combined  with  (3 -50).  Substituting  the  values 
of  q  and  n  into  these  equations^  we  obtain 


(2)    (k'|r|k)  =  ri  p  F(k)  { 


ih«i M/v-i   2k.(k-k')  -  (k-k')^ 


k-k' 


Ik-k'l^ 


.jfc  .       -k  -A  j» 
Let  us  now  compute  S^tk^P   and  g^Ck^+k^' k^) ,  We  recall  that 

where  ^(r)  is  the  wave  function  of  the  bound  electron  and  the  subscript  m 

indicates  the  state.  For  the  final  state,  we  are  using  a  plane  wave.  Hence 

,     j^    ~  ±kX'  T  +  ik^-r 
?  1  dr  e 
(Srt)' 


(M 


=  b{k^  -   kp, 


The  ground- state  wave  franc t ion  ^(r)   is 

Ut)   =  — ^T72  6""^^°   . 
Hence 

(5)    go^^i -^  %  -  V  =  iT-^f-jrv-F^   • 

rtjl+a^ik£+k^'-  k^l  j 

.A 

Substituting  (k)   and  (5)  into  (l),  the  k"  -integration  can  be  carried  through, 

yielding 
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1  2  ion  1       «   [l+a2j?.+?.-  ?j2j  2 

Cooiblnlng  (2)  and  (6)  and  substituting  into  (3.36)  we  obtain  for  the  differential 
scattering  cross  section 

(■^^   ^ion'^1^2  =  ^"P^ p/  ?---  -Tl 

ion  X  ^   -h  a  k-   sinh  it/a  2k-"k,'-k; 
o  1      L  '  o'  1  1  2'J 

6(E^+  Eg  -  E)dk^di^ 


1-*^  -^  \\x   ^    -*  -»•  i2r   21-^  -»■  -^  121 

IV4l  121^1-4-^21  toiV4-^2i  J 


The  total  cross  section  will  be 


<^t  =  I  /  qon  <^i^2 


'2' 
Before  we  proceed  to  the  evaluation  of  this  double  integral,  a  word  must  be 

said  about  the  limits  of  integration.  Due  to  the  5-function  appearing  in  the 

f 
integrand  of  (7)*  the  values  of  Ep  and  hence  k'  are  restricted  by  the  relation 

2 

2   -*     2 
k^  =  k^  -  1  -  k^''  . 

Hence  the  range  of  k'  will  have  to  extend  over  the  entire  space  of  k'.  We 

must  be  careful  in  dealing  with  the  &-function.  It  involves  the  square  of 

the  magnitude  of  k'  rather  than  k^.  This  is  dealt  with  by  using  the  relation 

SCx^-a^)  =  1^  (S(x-a)  +  5(x+a)). 

(This  point  is  dealt  with  in  Appendix  B.)  The  remaining  problem  is  the 
magnitudes  of  k ' .  Again  the  angular  integrations  can  be  deeilt  with  by 
substitution  of  variables.  This  is  done  in  Appendix  B.  The  range  of  k'  is 
fixed  by  conservation  of  energy.  Two  extreme  cases  occur  -  one  in  which  the 

hydrogen  atom  is  barely  ionized.  The  incident  electron  then  carries  off 

2    2 
the  entire  energy  and  k'  =  k^  -1.  At  the  other  end,  the  emerging  electron 

2 
carries  off  the  entire  energy  and  thus  k'  =  0. 

The  integration  is  carried  out  in  Appendix  B  and  leads  to 


2 


V 


^1 
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-15^1-^1  V  15^5 -^  3  ^^  +  "^^5- "^^6 


where  the  cr's  are  defined  in  the  appendix.  The  integrals  a'  cannot  be  carried 
out  analytically.  They  were  integrated  numerically.  The  resvilts  are  plotted 
in  Figure  1.  The  curves  axe  drawn  in  this  figure  and  are  labeled  according 
to  the  calculation  they  are  based  upon.  The  Born  approximation  was  computed 
numerically  by  Fite  who  used  the  Mott  and  Massey  formula  given  in  Section  2. 

The  maximum  occurs  at  20  volts  instead  of  98  volts  and  gives  a  fairly  large 

2 
value,  namely  1.2na  .  This  is  much  larger  than  the  experimental  value.  Our 

curve,   on  the  other  hand,   gives  a  maximum  of  0.55Jta     at  an  energy  of  86  volts 

and  is  in  fair  agreement  with  the  experimental  result. 

V/e  examine  now  the  behavior  of  the  cross  section  for  high  energies.  As 

iL,  ->  CO,  the  factor  1  +  (k  -  K)^  behaves  as  k.;  the  factor  e""'  (sinh  ^)— >  1 

and  thus  the  cross  section  behaves  a  l/k'  .   It  is  seen  that  it  falls  off 

very  rapidly  as  k^->  cD  . 

We  now  examine  the  behavior  at  threshold.   To  this  end,  we  go  back  to 

p 
equation  (7).  We  are  concerned  with  the  behavior  of  a  as  k^  =  1  +€  (6->0). 

Since  in  this  limit  k'  and  k'  are  both  small,  we  can  ^/rite, 

2Ve"''^^^  6(k'2+k'2-0 

0-.  ^dk'dk-  ^  ^-^ .^.^    r    ,^  I   ^  ,pU    ^{^2 

^°^    ^    2  4k5  sinh«/2kjL   |k^-k||[l+lk^-ki-k^|2p        ^    2 

We  carry  out  the  k'  integration;  this  leads  to 

Ihe  integral  over  k'  can  now  be  carried  throu^.  Ihe  limits  on  kJ  are  0 

1/2 
and  €  '  .  Integrating  leads  to 

01  =  const  € . 
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WaonlerL-'-  -I  derives  the  threshold  law  for  single  ionization  and  obtains 

_  1-123 
°t  **^ 

Massey  points  out  that  on  the  basis  of  the  available  evidence 

o:^  oce  . 

Our  resvilt  is  in  agreement  with  Massey 'sj  however,  it  must  be  remembered  that 
at  threshold  our  approximation  cannot  be  very  accurate. 

5.   1S-2S  cross  section 

A.  Direct  scattered  amplitude 

In  the  case  of  excitation,  the  R-matrix  takes  the  form 

(1)   (k'3ilR^|kO)  =  e^  I  dkg  ^(k2+k-k')g^(k^)(k^|r|k2) 

.A         -3k 

where  k  and  k'  are  the  initial  final  wave  niimbers  of  the  incident  electron 
respectively  and  n  and  0  are  symbols  representing,  respectively,  the  final 
and  initial  states  of  the  bound  electron  and 
(2a)  k^  =  2k •  -  k  -  k^^ 


(2b)         k^  =  k  -  k 


2 


Referring  to  Appendix  A,  we  have,  upon  the  use  of  (2), 

(3)     (4lHki)=-^^. 
^         ^  k-k'  l'^  ' 


(k-k')  -(k-k').(k-k2) 


In  oixler  to  carry  our  work  further,  we  must  make  one  more  approximation.   It 
is  clear  that  if  we  substitute  (5)  into  (l),  it  would  be  impossible  to  carry 
out  the  kp  integration.  Hence,  we  must  attempt  to  take  the  r-matrix  from  xmder 
the  integral  sign.  Chew L J  points  out  that  this  is  a  slowly  varying  matrix. 
We  will  therefore  substitute  for  k_  in  the  r-matrix  k  -  k ' .  This  approximation 
is  not  necessary  in  the  ionization  case.  The   reason  is  that  the  final  integral 
over  k^  involves  a  6-fvmction  of  kp  since  the  final  state  of  the  emitted 
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electron  is  represented  by  a  plane  wave.  If  we  had  used  exact  wave  fimctions, 
the  same  approximation  wovild  have  been  made.  Tlius, equation  (l)  beccanes 


w 


i/2a  k' 


it-  niHjfo) ,  -  i^ .  'g-^yy '      °  |a.vA)«;(-v-) 


We  recall  that 
(5a> 


g^(k2)  =  (2«) 


-5/2 


dr  e      ^(r) 


and 
(5b) 

Hence 


g^Ck+k^-k')  =  (2rt) 


■5/2 


^  -i(k+k  -k').r    ^ 
dr  e  ^(r) 


(6)    [dk2g^(k2)^(k+k2-k')  =  [d?fd?'J^(r)^(r')e^^^-'^-^'*^'^(2n)-^fdi2e  ^'   '""'' 


=  |d?0»(r)e^^^-^')-^(?). 


This  integral  is  evaliiated  in  Appendix  C;  one  then  obtains  for  the  direct 
13-23   cross  section 


(7)  (t',200|R^|k,100)  =  -  f  ■-  TJ^%.l   . 

[a-f4a^|k-k'rj^ 


(k-k')»(k+2k') 

(k-k-r 


i/2a  k' 
'  o 


B.  Exchange  scattered  ancplitude 

In  this  case  we  evalxiate  equation  (2.42)  which  is 


y 


(8)  (k'200|R^|k,100)  =  -  ie'='Jdk2gJk2)^QQ(k+k2-k')l 

where  ^ 

/  ,^  ik •  -r 

(9)  I.JS^e  ^   tf^(?). 


Referring  to  Appendix  A,  we  have 


I  = 


iffliN(k') 
2k'+k 


(k-k').(k+2k') 


(k+2k') 


.^2 


i/2a  k' 


Ihus 
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(kS200|R^^|k,100) 


4rte^N(k ' ) 
(k+^'  )^ 


(k-k ' ) .  (kH-2k ' ) 


"l/2a  k' 


(k+2k')' 


•^220^^2  ^^^"^'■'"V^^ 


or 


2^9/2«e\^(k')|k-k'|^ 


(10)  (k',200|R   |k,100)  = 

^         |2fi'+S|^|9+4a^|k-^'|2p 

C.  Scattering  cross  sections; 

The  scattering  cross  section  Is  given  by 


(k-g')-(g+2g') 
(g+2k')^ 


l/2a  k' 
'  o 


(11)    -=(^)^|^  ) 


\\t*if*  l\t-sf 


where  f  and  g  are  the  direct  and  exchange  aniplltiides .  This  can  also  he  written 
In  the  form 


(12)  cr    =  (^)2     |:  i 


|f|^+|g|2-  Refg* 


Substituting  the  values  of  f  and  g  from  (7)  and  (lO)  leads  to 


^■"■^^      °icX)->200 


kk'slnh 


-  [9+J+  |k-k '  I  ^Y     I        (k+^* )  (k+^'  )^ 


cos 


2k 


.loff^TT 


-*•    -*•  i2i 

k+2k'l 


k-k 


m2 


The  total  cross  section  Is 

(^^^  ^100->200  =  2«  I  ^100->200  ^^  ^®- 


nils  integral  was  ccnrputed  ninnerically  and  led  to  the  curve  given  in  Figure  2. 
Coaaparlson  with  the  experimental  results  of  Llchten  L-^-i  shows  that  the 
maximum  we  obtain  is  of  the  same  order  of  magnitude  but  occurs  at  a  higher 
energy.  Lichten's  maxljmim  occurs  at  11.7  e. v. 
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6.    1S-2P  transition 
A.  Direct  scattering 

The  main  difference  between  the  calciilation  for  the  2P  state  and  the  2S 
state  lies  in  the  evaluation  of  the  integral 


(1)    1=1  di2g^(k2)^(k+k2-k') 


since  the  final  state  is  not  the  same.  Hiis  integral  leads  to 


(2) 


^15/2  |e>_^, 


=  if^^ 


k-k' 


^a^|k-k'|^ 


Ihis  leads  to 


(3)        (kS210lRJk,100)  =  - 


2^'^/^3iJte^a  N(k') 


|k-k'|r9+^ 


a^ij^fp' 


(k-k')-(k+2k') 
(k-k')2 


i/2a  k' 
'  o 


1 


B.  Exchange  scattered  amplitude 

Again,  the  difference  between  the  2S  and  2P  lies  in  the  momentum  integral  (l). 


Hence  the  exchange  scattered  amplitude  is 

2^''/2«e\N(k • )  |k-k •  I    r  .,^,  ^    p 

(k)        (K',210|R   |k,100)  =  ^ 5 -^rr       ^^^  )'{y^+2k   ) 

^^         lk+2?'|^|9+^a2|^.^'|2  3  ^|^^^.|2 

C.  Scattering  cross  section 

The  scattering  cross  section  is  obtained  from  equation  (5-12) 
combined  with  (3)  and  (4).  This  yields 


Ti/2a  k' 


^^^        ^00->200 


-.  _16     -rt/k' 

9x2     ite     ' 


m  •  \t-P  l^sir^V  [^  l«'  1^]^  1      '"^  '■"'* 


(2k'+k)' 


cos 


/  i_  log  ikt£:£\ 


Ite  total  cross  section  would  then  be 
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It 
(6)   L  ioo->210  =  ^TT  j   ^ioO->210  ^^  ®*® 


Hiis  integral  was  computed  numerically  and  led  to  the  c\irve  given  in 

Figure  5.  Also  included  in  this  figure  are  the  first  Born  approximation 

calculation l- J,  the  second  Bom  approximation'--'  and  Flte's  experimental  result*-  -^  • 

Although  the  maximum  in  oxir  curve  is  shifted  toward  the  higher  energies, 

it  seems  to  us  that  it  is  a  much  better  result  than  either  of  the  approximations 

and  simpler  to  calculate. 

7.    Discussion 

Hie  cross  section  for  ionization  and  excitation  to  the  2S  and  2P  states 
of  hydrogen  have  been  calctilated  on  the  basis  of  the  'Impulse  Approximation." 
We  note  in  Section  5  that,  for  the  2S-state,  the  agreement  is  not  good 
although  the  magnitude  is  the  same.  For  ionization  and  excitation  to  the 
2P-state,  the  agreement  is  excellent  and  by  far  better  than  any  other  approximation 
method  that  is  available.   It  therefore  seems  that  it  is  possible  to  use  the 
Impulse  Approximation  for  calcxilating  cross  sections  in  the  scattering  of  three 
or  more  body  problems. 
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Appendix  A 

¥e  evaluate  here  the  integral  given  in  equation  (5.47) 


(1)    ^  =  jr  ^ F(-iD;l;ikr  -  i^.t) 


We  use  the  following  integral  representation  for  the  confluent  hyper geometric 

function  ^^^ 

(Of) 
(2)    F(-in:l;lkr  -  ik.?)  =  ^  J    u-^-^(u-l)^"e^^^  "  ^'^du  . 

1 

Substitution  of  (2)  into  (l)  yields 

(Of) 
(5)     1  =  2^  j    d^  ^"''"'^u-D^^Ku) 

1 

where 

{h)         l(u)  =   I  — -  expjiq-r  +  ikru-  ik-rul    . 

I(u)  is  a  standard  integral  that  can  be  carried  throvigh  to  yield 

q  -2q-ku 

(5) 


2rt 


(q.k)(u-u^) 


vlth        ^o  =  ^  /2q'k     . 

Thus 

(Of) 


(6)        1=   —r    f       u-^-^(u-i)^(u-u J-^du 


q-k        ^ 
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The  contotir  for  the  representation  of 
the  hypergeometric  function  is  given 
In  Figure  1.  Examination  of  the  integrand 
in  (6)  shows  that  we  have  branch  points 
at  0  and  1  and  a  simple  pole  at  u  .  Thus 
we  must  deform  the  contour  in  such  a  way 
as  to  be  able  to  include  the  pole  u  and 
not  the  branch  points.  Ttils   is  accom- 
plished by  the  contour  given  in  Figure 
2^     K     This  contour  is  now  clockwise. 
To  compensate  for  the  negative  sign 
thus  introduced,  we  put  a  negative  sign 
in  front  of  the  integrand.  Then  applying 
Cauchy's  theorem  we  get: 

«    -1   u  -1  . 
T  2«   1  /  o  xin 

q-h   o    o 


and  since 

we  get 
(7) 


%  =  -^^' 


2q.h 


I  = 


kn    rq^-2q-h  T 
n   L.   o      J 


FIG. 


^m 


FIG. 2 


which  is  the  integral  given  in  the  text. 
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Appendix  B 

We  evaluate  the  integral  occ^l^ing  in  the  total  ionization  cross  section, 
Dividing  01  as  given  "by  eqiiation  (U.8)  and  putting  h=m?:e=l,  we  have 

.83 


(1)   61  = 


2 
The  cross  section  is  expressed  now  in  units  of  3ta  ;  the  unit  of  length 

is  the  Bohr  radius,  and  the  unit  of  energy  twice  the  ionization  energy  of 

hydrogen. 


Let 


-V,      -!>.      -k. 

k'  =  k^  -  k^ 


Then: 

(2) 

Let 


,    ^i^[{  dkdk-  e'^/'^^^'l     y   (k^^k^^-k-^.H-2^.(k;k'))) 

^*  °  ^1  n         sinh  Tt/ik+r-i  k^k+k')2A+|k^.k-k'|^p 


(3)      k-k'  =  x[     [k  «  l/2(K+t) 
k+k'  =  Kl   ^  lk'=  l/2(K-X) 


then 


dkdk  •  =  ir  dXdK 
2^ 


Thus 


.    .     ^    5  jrx2-(i;t  .K-2k2-K2.2)ll 

(5)   Kic)  =:   dx   ^^'-^^.^    ^ -LL 

y         (K+X)^ 
The  angular  integration  can  he  carried  through  easily  leading  to 
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CD 


(6)    I(K)  =  kn   \     X' 


6{|(x2-[Uk^-K-2k^.K2.2])} 


■00 


(K^-X^)^ 


dX  . 


The  limits,  due  to  the  5-f unction^  will  run  over  the  entire  space. 
2 


Now  let 


a  =  iHt^-K  -  2kf  -  K^  -  2 


2 
^1 


Then 


8(J  (X^-a^))  =  U8(x2.a2) 


=  -(5(X+a)  +  5(X-a))  . 


Substituting  into  (6)  and  carrying  out  the  integration  yields: 

{    ^    ^      op      )  1/2 
h\k-k.  ■K-2kr-Yr-2\ 
(7)  i(K)  =  -5--i i ^ 


|l+|k^-K|2|2 


Then 


°t  = 


dVf^  e-"/^     ^^^i-^-^i-^:")"^^ 


2)T 


1     >         ^^^^^     )^      [l+1^4^-2k,Kx]^ 


^12 


dx 


Carrying  through  the  second  integral,  we  obtain 


£^'   .£S£^„.i^.,.MV^.i^ 


'*-Sf"^*;4^'"^*Sf '^"^'l""'""^"' 


where 


^1  = 


[l+(k^-K)2]5  ■  [l+(k^-K)2]5 


dK  e-^/^ 
'  K     sinh  rt/K 


°'2 


=  1 


/-Uk^K+K^+2k^+2 


V^ 


ifk  K+2k^+K^+2 

[i-KvK)T 


sinh  it/K  ^ 
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^5  = 


[l+(k^-K)2] 


■Jt/K 


sinh  Jt/K     j^5 


■n/K 


dK 

212        >  sinh  n/K  ^7 

A. 


/         /-4k^K+K^+2k^+2  Ak^K+2k^' 


2 
+k^+2 


[l+(k^-K)2] 


e-'^/^         dK 
'  sinh  jt/K  j^9 


o>  = 


tanh 


^  K  /+4k^K+2k^+K^+2                          -l^J  -^^iK+2k^+K^+2 
— ^ — — — r- +      -    tanh        — * ;:; — 'X 


2+2k?+K^ 


2+2k^+K^ 


e-^/^  ^  dK 
sinh  n/K  j^O 


The  limits  of  integration  are 


K  .  =  2k 
mln 


K    =  2k 
max 


1-Ar^ 


This  follows  from  the  discussion  in  Section  h. 
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Appendix  C 
We  evaliiate  here  the  integral  appearing  in  Section  5,  namely  eqxiation  6. 

1)   I„  =  |d?^^(?)e'(f-^''-^^(?).. 


( 

Now 

(2a)  |^(r)  =  (Wa5)-l/2  ^  ^^%   ^ 

For  final  states,  we  have 

(2b)  2S:  $200^r)  =  (2\a^)-l/2  (2-r/a^)e  ""     ^° 

-r/2a 


(2c)  2P:  1,,o(^  =   P-o]'"^"     r- 


1.  2S- state 


<"     ^--afTfe   ]-=''"''•' (^-/"o)' 


COS  e. 


■5r/2a, 


Let 


X  =  —  and  K  =  a  (k-k') 

a  o 

o 


then 


^200  "  rf7 


00 


2dxe-5^/2  (2-x)    (  e^  '°'  ®- 


i 


sin  ede  I     d0 
o 


CD 


{^) 


1  ^       ^        /o     \   -3x/2,    iKx     -iKx» 

I —     I       dx  x(2-x)e  -^  '    (e       -e         ) 


T^ix 


^^ 


o 
/ 


liX^-  X*)  -  (Xg  -  X*) 


y 


where 


(5a) 


(5b) 


oo 

o 
oo 


xdxe 


•(5/2  -   iK)x 


Xg  »   I     x2dxe-(5/2  -  iK)x 
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These  are  standard  integrals  and  yield 
(6a)    X  =  —^ ^ 

(6b)    X^  =   ^^ 


^  '  (5-2iK)^ 


Hence 


3-2^iK 


_  2  iK(g7- 
2   '^^     (9+*^K^)' 


(7b)   X,  -  X*  =  iMnd^  , 


Substituting  into  {h),   ve  get 

^15/2^      2^5/2,- ^.,2 


(8)    I 


{9-^Y^y    ~      (944jk-k'|^j 


2.  2P- state 


(5)        ^210  =  TiW^  I  "^  ^ 

2  '  «a^  j 


+i(k.k')-?  "^''/^^ 

^    '   e       cos  e 


The  same  substitution  is  carried  out  for  r  and  k-k'  as  for  Ip^y-,^  yielding 

oo  n 

1         (       3j     -5x/2     /     iKx  cos  e  _     .     ^_ 

210  ~  "5/2      I     ^  *^®  J   ®  °°^  ®  ^^'^  ^®   • 


o  o 


These  are  simple  integrals  and  can  be  carried  th^o^^gh  to  yield 

(10)     I     -  ?ig"'^^|g-g'l 
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